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The origin of the lead-rich stars in Galactic halo: investigation of the 
model parameters for the s-process 
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ABSTRACT 

Several stars at the low-metallicity extreme of the Galactic halo show large spreads of [Pb/hs]. 
Theoretically, a s-process pattern should be obtained from an AGB star with fixed metallicity 
and initial mass. For the third dredge-up and the s -process model, several important properties 
depend primarily on the core mass of AGB stars. Ziilstra ( 2004) reported that the initial-final- 
mass relation steepens at low metallicity, due to low mass-loss efficiency. This perhaps affects 
the model parameters of the AGB stars, e.g. the overlap factor and the neutron irradiation time, 
in particular at low metallicity. The calculated results show indeed that the overlap factor 
and the neutron irradiation time are significantly small at low metallicities, especially for 
3.OM AGB stars. The scatter of [Pb/hs] found in low metallicities can therefore be explained 
naturally when varying the initial mass of the low-mass AGB stars. 
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1 INTRODUCTION 

The elements heavier than the iron peak are made through 
neutron capture via two principal processes: the r-process (for 
rapid p rocess) and the s- process (for slow neutron capture 
process)! Burbid ge et aljfl957i) . The observations have confirmed 
that, indeed, Asymptotic Giant Branch (AGB) stars show over- 
abundances of elements heavier than iron at their surface 
JSmith & Lambert 1990), which clearly indicate that the s-process 
takes place during the AGB phase in the evolution of low- and 
intermediate-mass stars (0.8^M(Mq)^8). Low-mass AGB stars 
are usually thoug ht as the main nucl ear productio n site of the s- 
process elements iGallino et all 1 9981: iLugaro et alJl2003tlHerwid 
2004). By now, the generally favoured s-process model is as- 
sociated with the partial mixing of protons (PMP) into t he ra- 
diative C-rich layer s during thermal pulses jStra niero et al1ll995t 
lGallinoetal1ll998l l2003t Istraniero et alj|20oj) . PMP activates 
the chain of reactions 12 C(p,7) 13 N(/3) 13 C(a,n) 16 which likely 
occurs in a narrow mass region of the He intershell (i.e. 13 C- 
pocket) during the interpulse phases of an AGB star. The s-elements 
thus produced in the deep interior by successive neutron captures 
are subsequently brought to the surface by the third dredge-up. 
Using the primary-like neutron source ( 13 C(a,n) ls O) and start- 
ing with a very low initial metallicity, most iron seeds are con- 
verted into 208 Pb. When the third dredge-up episodes mix the 
neutron capture products into the envelope, the star will appear 
s-enhanced and lead-rich. If the standard PMP scenario holds, 



all s-process-enriched AGB stars with metallicities [Fe/H] 1 ^ -1.3 
are thus predicted to be lead(Pb) stars ([Pb/hs] ^ 1, where hs de- 
notes the 'heavy' s-process elements such as Ba, La, Ce), inde- 
pendently of their initial mass and metallicity {Gallino et aljl998t 
iGorielv & MowlavifcOOOLlGorielv & Siesj200ll) .~ 

The first three such lead stars (H P 187861, HD224959 , 
HD196944 ), hav e been later confirm ed bv IVan Eck et alj feOOll) . 
At the same time. I Aoki et al J i20 1 ) found that the slightly more 
metal-deficient stars LP 625-44 and LP 706-7 are enriched in 
s-elements, but cannot be considered as lead stars ([Pb/Ce]< 
0.4), in disagreement with the standard PMP predictions. Re- 
cently, more spectroscop i c data of s-ri ch and lea d-rich stars are re- 
ported lAoki et alfeOoilCohen et alJ|2003l iLucatello et alj|2003l 
[ Van Eck et alj|2003t Ijohnson & Boltell2002L 12004 ISivarani etal] 
2004 ). The large observation data spreads of [Pb/hs] are strong indi- 
cation to suspect a large intrinsic spread of integrated neutron irra- 
diations. In order to explain the spreads of [Pb/hs], a large spread of 
13 C-p ocket efficiencies is subsequently proposed by Straniero et al. 
(2005). However, it should be stressed here that the predict ions of 
the st andard PMP scenario are rather robust iGorielv & Mowlavll 
l200Ct) . In the framework of the PMP scenario, there is no obvious 
degree of freedom that could be used to reduce the lead production 
in low-metallicity AGB stars (Van Eck et al. 2003). At present, the 
physical explanation for the different 13 C-pocket strengths, which 
perhaps should not be consistent with the primary na ture of the 
neutron source, is not yet found iRevniers et alj|2004l) . Thus the 
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1 where [X/Y]=log(N x / Ny) - log(N x / Ny ) Q , refers to the Solar Sys- 
tem abundances. 
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fundamental problems, such as the formation and the consistency 
of the 13 C-pocket, the neutron exposure signature in the interpulse 
period, currently exist in the models of AGB stars. 

For the third dredge-up and the thermal pulse model, sev- 
eral important properties depend primarily on the core mass M c , 
while the dependence on the other stellar parameters is negligible or 
margi nal jlberll9 77l:lGroenewegen & de Jong|l993tlKarakas et alj 
2002). Ziilstra 1 2004) reported that the initial-final-mass relation 
steepens at low metallicity, due to low mass-loss efficiency. This 
may cause the degenerate cores of low-Z, high-mass AGB stars to 
reach the Chan dresekhar mass, leading to an Iben & Renzini-type- 
1.5 supernova (Iben & Renzini 1983). On the other hand, this can 
obviously affect model parameters of the AGB stars, e.g. the over- 
lap factor r, which is the fraction of material that remains to ex- 
perience subsequent neutron exposures, and the neutron irradiation 
time At, in particular at low metallicity. 

In this paper we will present a calculation which indicates that 
at low metallicity, large core mass of AGB stars may allow the over- 
lap factor and the duration of neutron irradiation to reach small val- 
ues for a fixed initial mass of AGB stars. This will affect the charac- 
ters of the s-process nucleosynthesis, and can explain the observed 
abundance pattern of lead-rich stars. The next section discusses the 
model parameters of AGB stars. In section 3, we discuss the char- 
acteristics of the s-process at low metallicity and the possibility of 
lead stars. Finally, in section 4 we summarize the main conclusion 
that can be drawn from such an analysis. 



2 MODEL PARAMETERS OF AGB STARS 

There are four parameters in the parametric model o flHoward et alJ 
( 1986) on s-process nucleosynthesis. They are the neutron irradia- 
tion time At, the neutron number density N n , the temperature 7g(in 
units of 10 9 K) at the onset of the s-process, and the overlap factor 
r. Combining these quantities, we can obtain the neutron exposure, 
Ar=N„VTAt, where vt is the average thermal velocity of neu- 
trons at Ty. The temperature is fixed at a reasonable value for the 
13 C(a,n) 16 reaction, 7g=0. 1, for these studies. Using the initial- 
final m ass relations as function of metallicity presented bv lZiilstnl 
120041) . the effects on the parameters can be derived. 
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Figure 1. The overlap factor of different initial mass AGB stars, as function 
of metallicity. 

which is significantly small at low metallicities, especially for 3Mq 
AGB stars. In AGB stars with initial mass in the range M = 1.5 ~ 
3.0Mq, the core mass M c lies between 0.6 and 1.4Mq at [Fe/H]= 
-2.5. According to the equation ( 1), the corr esponding values of r 
will range between 0.8 and 0.13. iGallino et al J \ 1 9981) have found 
an overlap factor of r = 0.4 ~ 0.7 in their standard evolution model 
of low-mass AG B stars at solar m etallicity, which lies in our pre- 
dicted range of r. Aoki et al. ( 2001) have reported an overlap factor 
of r ~ 0.1, found for the best fit to metal-deficient AGB stars that 
produced the abundance patterns of LP 625-44 and LP 706-7. In an 
evolution model of AGB stars, a small r may be realized if the third 
dredge-up is deep enough for s-processed m aterial to be diluted by 
extensive admixture of unprocessed material. llCarakas et alj 1 2002) 
have found that the third dredge-up is more efficient for the AGB 
stars with larger core mass. Taking account of the core-mass depen- 
dence, the wide range of r-values of the lead enhanced stars can be 
explained naturally by the wide range of core-mass values of AGB 
stars at low metallicity. 



2.1 The overlap factor 

For the AGB model, the overlap factor r is a fundam ental parame- 
ter. An analytical formula was given by Iben ( 1977) as a function 
of the core mass M c in the range 0.6^M C ^ 1.36: 

r = 0.43 - 0.795(M C - 0.96) + 0.346(M C - 0.96) 2 . (1) 

We can obtain an initial-final mass relation as a function of metal- 
licity, by fitting results of Ziilstra (2004): 

M c = A(Z) + B(Z)M, (2) 
where 

A(Z) = 0.46449 + 0.03279log-^- + 0.000U(log-^-) 2 , (3) 

Zq Zq 

and 

B(Z) = 0.08729 - 0.04851/og^- + 0.00468(Zoy^-) 2 , (4) 

Zq Zq 

which is valid for -4^log-^^0. Combing the equations (1), (2), 
(3) and (4), we obtain the overlap factor as a function of the ini- 
tial mass and metallicity. The overlap factor is shown in Fig. 1, 



2.2 The neutron exposure 

IGallino et alJ 1 19981) have pointed out that the neutron density is 
relatively low, reaching ~10 7 cm -3 at solar metallicity. Since the 
13 C neutron source is of primary nature, the typical neutron den- 
sity in the nucleosynthesis zone scales roughly as 1/Z 6 , from Zq 
down to 1/50Z(7) . At lower metallicities, the effect of the primary 
poisons prevails (Gallino et al. 1999; Busso et al. 1999). 

There is a possibility for the synthesis of s-process elements in 
the AGB st ars, i.e., wit h nucleosynthesis taking place during ther- 
mal pulses <Aoki et alj200ll) . In this case, the neutron irradiation is 
derived primarily by the 13 C(a,n) 16 reaction, with a minor con- 
tribution from the marginal burning of 22 Ne. 

However, in the s-process scenario that invokes radiative 13 C- 
burning, the nucleosynthesis mostly occurs during the relatively 
long interpulse period, in a thin radiative layer at the top of the 
He intershell (i.e., the 13 C-pocket model). A second neutron burst 
giving rise to a small neutron exposure is released by the marginal 
activation of the 22 Ne neutron source in the convective thermal 
pulse. The neutron irradiation time A? should be close to the in- 
terpulse period at low metallicity due to the combination of two 
reasons. The first is that the higher neutron density can lead to 
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Figure 2. The neutron exposure of different initial mass AGB stars, as function of metallicity. 



longer neutron irradiation time, and the second is that the shorter 
interpulse period is expected for larger core-mass of AGB stars 
iGroeneweeen & de .Tonglll993l) . Therefor, adopting the interpulse 
period as the neutron irradiation time will have a smaller effect on 
the low-metallicity stars of interest here than that on stars of solar 
metallicity. 

In our calculation, the neutron irradiation time At is adopted 
respectively as follows: 

Case A: the duration of the thermal pulse, where the core- mass- 
duratio n of the convective shell relation is adopted from llbenl 
ll977l) . 

Case B: the interpulse perio d, where the core-mass-in terpulse pe- 
riod relation is adopted from Boothro vd & Sackmannl ll988l) . 
We choose respectively Ar=0.08mb -1 at [Fe/H]= -0.3 in 
case A, which corresponds to a mean neutron exposure 
r =0.296(r g /0.348) 1/2 mb~ 1 , and AT =02mb_^for 3.0M n AGB 
stars with solar metallicity in case B iGallino et al. 1998). Using 
the initial-final mass relations given by Ziilstra (2004) and the neu- 
tron irradiation time At, we can obtain the neutron exposure Ar 
as a function of metallicity and initial mass (see Fig. 2). The trend 
shown in Fig.2 (Case A and Case B) can be understood as follows: 
At is proportional to the neutron number density N n and the neu- 
tron irradiation time At, where N„ is expected to increase with de- 
clining metallicity. However, Ar declines with declining metallicity 
due to the increasing of stellar core mass, which directly leads to 
a decline of Ar at very low metallicity, especially for 3M0 AGB 
stars. 

Based on the primary nature of the 13 C neutron source, the 
value of Ar will reach about 6.3 mb -1 around [Fe/H l^ -2.5 for 
the case of radiative 13 C-burning iGallino etall ll 999). Our result 
of case B for the 1.5M© AGB stars is close to the above value. Be- 
cause the neutron irradiation time At is shorter for the larger AGB 
stars, the neutron exposure Ar should be smaller too. lAoki et all 
1200 ll) have reported a neutron exposure, Ar ~ 0.7 mb _1 for 
metal-deficient stars LP 625-44 and LP 706-7, which is in the range 
of our calculated results for the both cases. The results shown in 
Fig. 2 imply that the wide range of Ar can be obtained naturally 
by considering the dependence of the irradiation time on the core 
mass . Since the Pb abunda nce is very sensitiv e to the neutron ex- 
posure I Galli no et all2003tlLugaro et all2003l) . large variations of 
the [Pb/hs] ratio could be expected. 



3 DISCUSSION 



3.1 The Case A 

The first model for 13 C-burning in AGB stars assumed that the 
neutrons were released in convective conditions (Hollo well&Ibenl 
1988; Kappeler et al. 1990). In such calculations, a repeated neu- 
tron exposure was achieved thanks to partial overlapping of mate- 
rial cycled through several thermal pulses. The s-process mecha- 
nism could be approximated by an exponential distribution of neu- 
tron exposures oc exp(— r/ro), where the mean neutron exposure is 
given by ro= -Ar/lw. The final abundance distributions depend 
mainly upon ro. 

In order to investigate t he efficiency of the s -process site, 
[Pb/hs] is particularly useful tStran iero et alj|2005l) . There have 
been many theoretica l studies of s-process nu cleosynth esis in 
low-mass AGB stars iDelaude etafll2004 llwamoto et alJl2004 
IStraniero et al]|2005l) . Unfortunately, the precise mechanism for 
chemical mixing of protons from the hydrogen-rich envelop 
into the 12 C-rich la yer to form 13 C-pocket is still unknown 
iRevniers et aljl2004 . This makes it even harder to understand a 
large spread of [Pb/hs] found in carbon-rich, metal-deficient stars. 
It is an interesting exercise to investigate the effect of the pa- 
rameters presented above upon the s-process efficiency of AGB 
stars. For this purpose , we have used t he simple analytical formu- 
lation IClavton & Rassbachlll967tlClavton & Wan! 1974 without 
depending on any specific stellar model, with many of the neutron- 
capture rates updated I Bao et al. 2000), to study what physical con- 
ditions are possible to reproduce the observed abundance pattern 
found in the metal-poor stars. The variation of the logarithmic ratio 
[Pb/hs] with metallicity is shown in Fig. 3a, where solid lines repre- 
sent respectively results of different initial mass of AGB stars. As a 
comparison, spectroscopic measurements (filled squares) of C and 
s-rich metal-poor stars are reported. Because of the uncertainties 
relate d to the formation mechanism of the 13 C-pocket iBusso et alj 
2001), a large spread of 13 C-pocket efficiencies has been proposed 
bv lStraniero et alJl2005t) in order to explain the spre ads of [Pb/hs], 
which has proved to be effec t ive for several purp oses iGallino et alJ 
Il998l iTravaslio et alJll999t IBusso et alj|2001l) . The results (plot 
lines) predicted bv IStrantox^HuTT2005l) for their standard case 
(hereafter ST), STx 2 and ST/75 case a re also presented respec- 
tively. The ST case iGallino etallll 9981) was shown to reproduce 



4 Wenyuan Cui, and Bo Zhang 




Figure 3. (a) Predicted [Pb/hs] versus [Fe/H] for different initial stellar mass choices (solid lines). Plot lines represen t the results predicted bv lStraniero et alj 
i 2005 ) for their ST, ST X 2 and ST/75 case respectively. Spectroscopic data of s-rich and lead-rich stars adopted from Delaude et al. 1 2004 ) are also included 
for comparison, (b) The ratio of lead stars over the low mass AGB stars (1.5-3.0Mq), as function of metallicity. 




Figure 4. (a) Predicted [Pb/hs] versus [Fe/H] for different initial stellar mass choices (solid lines). Plot lines represen t the results pre dicted by Stranier o et alj 
1 2005) for their ST, STx2 and ST/75 case respectively. Spectroscopic data of s-rich and lead-rich stars adopted from Delaude et ai] <2004 are also included 
for comparison, (b) The ratio of lead stars over the low mass AGB stars (1.5-3.0Mq), as function of metallicity. 



remarkably well the abundance distribution of the solar main com- 
ponent at metallicity Z=Zq/2 i Arlandini et al. 1999). 

Assuming a Salpeter IMF with a= -2.35, the ratio of lead stars 
compared to the low mass AGB stars' (1.5-3.OM0) is shown in Fig. 
3b. 



3.2 The Case B 

We use the parametric approach based on th e model of low-mass 
AGB stars computed bv lGallino et alj fl998l) assuming that all the 
pulses are identical. In this case, the neutron exposure is not well 
approximated with the exponential distribution, and the final s- 
process abundance distributions depend mainly upon the neutron 
exposure At, the mass fraction of 13 C-pocket in the He intershell 
q (adopted as 0.05) and overlap factor r. The adopted initial abun- 
dances of seed nuclei lighter than the iron peak elements were taken 
to be the solar-system abundances, scaled to the value of [Fe/H]. 
Because the neutron-capture-element component of the interstel- 
lar gas that formed very metal-deficient stars is expected to consist 
of mostly pure r-process elements, for the other heavier nuclei we 



use th e r-process abundances of the solar system l Arlandi ni et alj 
1999), normalized to the value of [Fe/H]. We carry out s-process 
nucleosynthesis calcul ation by means of an extensive reaction net- 
work described earlier dLiang et alj200d) . The results are shown in 
Fig. 4. 

It results that for very metal-poor stars, large spreads of 
[Pb/hs] are predicted for both cases. The agreement of the results 
with the observations provides a strong support to the validity of the 
parameters adopted in this work. Here, the logarithmic ratio [Pb/hs] 
shows a complex trend versus [Fe/H], due to the conjunct effect of 
the overlap factor and neutron exposure. Although the parameters 
of the two cases are different in part, the ratios of lead stars obtained 
in this work are very similar. At [Fe/H] ~ -2.0, the ratios of both 
cases all reach the maximum value around 0.55. The results ob- 
tained in this paper are the evidence that maybe the well established 
theories of the s-process nucleosynthesis, as it works in solar-like 
metallicities, may not work well at extremely low metallicity be- 
cause of greatly changed parameters. In fact, such r esults co uld be 
mainly led to by the new initial-final mass relations l Ziilstra 2004), 
i.e. the inefficient mass loss of the AGB stars at low metallicity. 
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4 CONCLUSION 

Theoretically, a s-process pattern should be obtained from an AGB 
star with fixed metallicity and initial mass. Taking account of the 
core-mass dependence, the large intrinsic spread of integrated neu- 
tron irradiations for the AGB stars at low metallicities is obtained, 
then the scatter of [Pb/hs] such as found in low metallicities can 
therefore be explained naturally when varying the initial mass of 
the AGB stars. Based on the relation of the overlap factor with the 
initial mass of the AGB stars, we can speculate that the Pb stars are 
polluted by low mass AGB stars (e.g. 1.5-2.5M0) and the non-Pb 
stars ([Pb/hs]<l) are polluted by larger mass AGB stars. We re- 
mind the reader that though the parametric approach is still useful 
to interpret observation data, it does not refer to detailed stellar evo- 
lution models. Obviously, a more precise overlap factor-core mass 
law and a more precise neutron irradiation time-core mass relation 
still have to wait for new models of nucleosynthesis in AGB stars. 
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